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OBJECTIVES 

The  objective  stated  in  the  application  was  to  carry  out  simulation  studies  on  the  5 
component  Vitreloy-1  metallic  glass.  This  was  to  include  the  characterisation  of 
interaction  potentials  and  the  running  of  preliminary  simulation  calculations.  One 
specific  goal  was  to  compare  the  outcomes  of  these  calculations  with  the  prediction  of  a 
crystal-based  structural  approximation  to  the  amorphous  state  recently  proposed  by 
Miracle.  The  project  was  originally  designed  to  cover  the  developmental  stages  of  a  more 
extended  study  of  an  important  glass  forming  system.  Soon  after  the  start  of  the  grant,  it 
became  uncertain  whether  support  could  be  ensured  for  this  extended  plan.  Furthermore, 
we  realised  that  a  key  feature  of  the  project  -  the  examination  of  the  Miracle  structural 
approximation  to  the  amorphous  alloy  -  could  be  effectively  treated  within  the  period  of 
the  grant  using  a  variety  of  hard  sphere  packing  protocols.  As  a  result  of  these 
considerations,  we  modified  the  objectives  to  the  following. 

The  objectives  of  the  work  reported  below  was  to  a)  explore  the  properties  and  stability 
of  the  crystalline  approximat  for  amorphous  alloys  proposed  by  Dr  Miracle  using  binary 
and  ternary  hard  sphere  mixtures,  b)  to  clarify  the  compositional  bounds  of  chemical 
order  in  amorphous  packings,  and  c)  to  establish  the  utility  of  hard  sphere  packings  as 
structural  representations  of  more  realistic  models  of  metallic  glasses. 


ACTIVITIES 


The  following  activities  have  taken  place  as  part  of  this  project. 
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1)  Prof.  Harrowell  and  Dr.  Fernandez  visited  the  Wright  Laboratories  in  Dayton,  Ohio 
during  Feb  7-11,  2005.  This  visit  included  extensive  discussions  with  Dr  Miracle  on  the 
proposed  research,  a  seminar  presented  by  Prof.  Harrowell  on  recent  research  in 
modeling  dynamics  in  amorphous  alloys,  discussions  with  members  of  Dr  Miracle’s 
group  and  a  visit  to  Dr  Miracle’s  lab. 

2)  Funds  from  the  grant  were  used  to  employ  two  research  assistants,  A.N.Alcaraz  and 
R.S.Duhau  from  CNEA  in  Buenos  Aires  from  Aprill,  2005  to  January  31,  2006.  These 
workers  performed  the  bulk  of  the  calculations  described  below  under  the  supervision  of 
Dr.  Fernandez. 

3)  On  Dec.  2,  2004,  Dr  Miracle  visited  the  University  of  Sydney  where  he  presented  a 
talk  on  modeling  the  structure  of  metallic  glasses  and  discussed  the  project  with  Prof. 
Harrowell. 


OUTCOMES 

1)  Computational  characterisation  of  the  Miracle  periodic  approximate  (MPA)  for 
binary  and  ternary  mixtures  of  hard  spheres. 


1.1)  Algorithm  Development.  The  algorithm  used  is  a  modification  of  one  described 
by  Clarke  and  Wiley  [Phys.Rev.B  35,  7350  (1987)]  which  efficiently  generates  the 
closest  packing  configuration  in  a  system  of  hard  spheres.  The  innovation  in  this  project 
was  to  use  as  starting  configurations  an  array  of  solute  particles  (labelled  B,  and  C,  in  the 
case  of  the  ternary  mixture)  on  the  sites  of  a  crystalline  FCC  lattice  (B  particles)  and  in 
the  octahedral  interstitial  sites  (C  particles)  in  accord  with  the  Miracle  periodic 
approximate  (MPA).  The  lattice  was  then  expanded  to  accommodate  the  random 
insertion  of  the  solvent  (A  particles).  This  procedure  ensured  that  solute-solute  contact 
was  avoided  initially  and  thus  introduced  chemical  order  into  the  hard  sphere  system  for 
which  there  is  no  such  ‘natural’  inclination.  The  relaxation  of  the  resulting  configuration 
resulted  in  the  configurations  that  we  refer  to  here  as  “constrained”.  “Unconstrained” 
configurations  where  generated  by  placing  particles  in  random  positions  in  the  simulation 
cell  and  maximising  the  density.  The  packing  fraction  r)  of  a  system  with  n  components 
and  no  overlap  (5=0)  is  defined  as: 


p  = 


47T 


i= 1 


V 


3 

where  V =nxnxnxa  is  the  volume  of  the  simulation  cell. 


As  a  validation  of  the  algorithm,  a  system  of  hard  spheres  of  equal  size  was  studied. 
The  maximum  packing  fraction  obtained  by  our  algorithm  was  0.641,  in  good  agreement 
with  the  accepted  value  for  random  close  packed  spheres  of  0.64. 
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1.2)  Properties  of  the  MPA.  Two  binary  (f?B/f?A=0.80  and  1,12)  and  one  ternary 
(f?B/f?A=0.80  and  RJ  RA={)  .62)  systems  with  various  compositions  were  studied.  The 
radius  ratio  chosen  correspond  to  definite  theoretical  coordination  numbers:  f?B/f?A=0.62 
corresponds  to  ZBA=8,  f?B/f?A=0.80  to  ZBA=10,  and  i?B/i?A=1.12  to  ZBA=15,  where  Z,,  is  the 
number  of  spheres  of  type  j  around  a  sphere  of  type  i. 

The  configurations  chosen  for  these  studies  were  those  with  the  highest  r)  values 
obtained  during  each  relaxation  process.  Pair  correlation  functions  gi/r)  and  coordination 
numbers  Z,y  were  calculated.  The  number  of  i-j  bonds  in  any  system  can  be  written  as 
NjZij=NjZjj,  then  it  is  possible  to  write  Zj^Z^XJX,,  where  Xj=NlIZNl  is  the  fraction  of  the 
component  i.  Therefore,  there  will  be  only  three  independent  coordination  numbers  (ZAA, 
ZBA  and  ZBB)  for  a  binary  system  and  six  (ZAA,  ZBA,  ZBB,  Zqa,  Zqb  and  Zee)  for  a  ternary. 

For  calculating  Zy  it  is  necessary  to  define  the  maximum  distance  r,j  for  which  two 
spheres  i  and  j  are  considered  nearest  neighbours.  This  distance  is  defined  as  the  one  at 
which  the  function  gi/r)  has  its  first  minimum.  The  r,y  distances  found  are  shown  in  Table 
I.  They  depend  on  the  size  ratio  but  are  independent  on  the  system  composition. 


binary  system 


Rb/Ra 

r aa 

>'AB 

7'BB 

0.8 

2.7 

2.5 

2.25 

1.12 

2.7 

2.9 

2. 7/3. 0’ 

ternary  sistem 


Rb/Ra:R(;/Ra 

r AA 

tab 

I'AC 

>'BB 

>'BC 

rcc 

0.8;0.62 

2.7 

2.5 

2.3 

2.2 

2.0 

1.7 

*value  corresponding  to  the  unconstrained  case 

Table  I:  First  minimum  distances  ri;-  (in  units  of  RA)  for  the  functions  gi/r)  in  the  studied  cases. 

The  sizes  of  the  systems  studied  ranged  from  1.0x10  to  6.8xl0~  spheres.  Most  of  them 
have  4x5x5x5=500  spheres  of  type  B  and  a  variable  number  of  spheres  of  type  A, 
according  to  the  value  of  the  chosen  fee  lattice  parameter  a.  The  composition  range  was 
Xb~0.07-0.50. 

Crystalline  order  in  the  solvents  was  found  around  two  composition  values,  0.125  and 
0.25.  The  crystal  structures  were  identified  as  Ca7Ge  (XB=0.125)  and  AlFe3  (XB=0.25,  3). 
The  presence  of  these  fully  ordered  structures  clearly  represents  a  complication  in  the  use 
of  the  MCA  in  describing  amorphous  states. 

The  packing  fraction  p  for  both  systems  with  either  size  ratio  is  plotted  as  function  of 
Xb  in  Figure  1.  Except  for  the  crystal  compositions,  there  are  no  important  differences 
between  constrained  and  unconstrained  systems.  Figure  2  shows  the  coordination 
numbers  Z(/  as  a  function  of  XB  for  both  systems.  The  peaks  and  cusps  at  AB=0. 1 25  and 
0.25  correspond  to  the  formation  of  crystal  structures.  It  can  be  seen  that  the  average 
coordination  numbers  ZAA  and  ZBA  are  similar  to  the  corresponding  theoretical 
coordination  numbers  in  the  low  composition  region  (ZAA«12,  ZBA«10  for  Rb/Ra= 0.8  and 
ZBA«  14-15  for  Rb/R\=\  .  1 2).  It  is  interesting  to  mention  that  cluster  sizes  distributions 
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were  rather  narrow  in  most  of  the  systems.  For  example,  around  50%  of  the  B  spheres 
have  Zba=10,  while  20%  have  Zba=9  and  another  20%  have  Zba=1  1  for  the  constrained 
systems  with  f?B/f?A=0.8  atXB<0.30. 


Figure  1:  Packing  fractions  rj  for  the  binary  systems  studied:  a)  i?B/^A=0.80,  b)  _/?b/^a=1- 12.  The  full  line 
and  symbols  (dashed  line)  correspond  to  the  constrained  (unconstrained)  system. 


Figure  2:  Coordination  numbers  Ztj  for  the  binary  systems  studied:  a)  7?B/7?A=0.80,  b)  R^/Ra=  1.12.  The  full 
line  and  symbols  (dashed  line)  correspond  to  the  constrained  (unconstrained)  system. 

1.3)  Stability  of  the  MPA.  Although  in  its  original  form  the  model  considers  that  the 
distortions  involved  destroy  the  long  range  order,  the  structure  analysis  is  much  easier 
when  the  B  spheres  are  assumed  to  occupy  fee  lattice  positions.  To  quantify  the 
magnitude  of  the  distortions,  the  displacement  suffered  by  the  B  spheres  from  the  initial 
perfect  fee  lattice  were  calculated.  Figure  3  shows  the  square  displacement  distributions 
for  the  whole  range  of  compositions  for  both  binary  systems  studied  and  Figure  4  the 
corresponding  mean  square  displacement.  Again,  the  cusps  and  peaks  are  due  to  the 
formation  of  crystal  structures.  In  general,  it  can  be  seen  that  the  displacements  from  the 
ideal  fee  lattice  are  of  the  order  of  f?A  or  smaller,  except  for  compositions  larger  than 
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around  XB=0. 3  where  the  number  of  B-B  contacts  increases  rapidly  and  the  subjacent  fee 
lattice  completely  falls  apart. 


1(T  10‘5  1(T  10'J  1(T  10'1  10° 

squared  displacement 


(a) 


io-5  io4  io-3  io'2  io''  io°  io1 
squared  displacement 


(b) 


Figure  3:  Squared  displacement  distribution  of  B  spheres  with  respect  to  their  ideal  fee  positions  for  the 
binary  hard  spheres  systems  studied  as  a  function  of  the  composition  XB:  a)  RB/RA=0.8,  b)  Rn/RA^\ .  12.  The 
squared  displacements  are  given  in  units  of  the  fee  lattice  parameter  a. 


(a)  (b) 

Figure  4:  Mean  squared  displacements  (MSD)  for  the  binary  hard  spheres  systems  studied  as  a  function  of 
the  composition  XB :  a)  RB/RA= 0.8,  b)  7?B/7?A=1.12.  The  MSD  is  given  in  units  of  the  radius  of  the  fee  lattice 
parameter  a. 

A  number  of  valuable  observations  relating  to  these  results  can  be  made. 


i)  We  find  that  the  original  solute  ordering  is  retained  upon  density  maximisation  up  to 
xb  =  0.19  for  R  =  0.8  and  R  =  1.12  (see  Figs  3  and  4).  Over  the  same  composition  range 
these  systems  have  a  density  within  1%  of  the  random  packing  case  and  of  configuration 
numbers  Zaa  very  close  to  that  observed  in  the  random  packing  case.  These  results 
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support  the  basic  contention  of  the  MCA  that  one  can  produce  a  very  reasonable  packing 
result  while  retaining  a  crystalline  organisation  of  the  solutes. 

ii)  The  presence  of  crystalline  phases  is  a  difficulty  for  the  MCA  in  that  they  create 
compositional  ‘blindspots’  where  we  cannot  extract  any  useful  information  about 
amorphous  packing.  Here  the  use  of  a  random  close  packed  structure  rather  than  fee  for 
the  alpha  solutes  would  be  of  benefit. 

1.4)  The  MPA  for  a  Ternary  Alloy.  The  total  number  of  spheres  in  these  systems 
range  from  1.7xl03  to  8.7xl03.  The  number  of  B  and  C  spheres  were  the  same  (i.e., 
Xb=Xq),  500  of  each  type.  Therefore,  the  composition  range  was  Ab~0.06-0.25. 

As  in  the  binary  system  case,  there  are  two  compositions  for  which  crystal  formation 
was  found.  ForXB=0.25,  the  structure  formed  corresponds  to  AlMnCu2  (L2i):  A1  and  Mn 


Figure  5:  Packing  fraction  r|  for  the  ternary  system  studied,  _/?b/^a=0.80  and  Rc/Ra=0.62.  The  full  line  and 
symbols  (dashed  line)  correspond  to  the  constrained  (unconstrained)  system. 


occupy  NaCl  positions  while  Cu  sits  on  tetrahedral  interstices  of  the  fee  lattice.  The 
AlFe3  can  be  obtained  out  of  this  structure  if  the  Mn  atoms  are  replaced  by  Cu. 

For  compositions  close  to  XB=0.125,  the  crystal  formed  could  not  be  identified  with 
any  known  structure.  The  packing  fraction  p  as  a  function  of  the  composition  X\>=Xq  is 
shown  in  Figure  5.  It  can  be  seen  that,  except  for  the  crystal  compositions,  no  essential 
difference  exist  between  the  values  obtained  for  constrained  and  unconstrained  systems. 

The  coordination  numbers  Z(/  are  shown  in  Figure  6.  ZA a,  Zba  and  ZcA  tend  towards 
their  theoretical  values  12,  10  and  8,  respectively,  when  the  composition  decreases.  The 
behaviour  of  ZBB,  ZcB  and  Zee  with  XB  is  related  to  the  NaCl  lattice  relaxation  in  the 
environment  of  A  spheres.  In  the  NaCl  lattice,  the  distance  between  spheres  of  the  same 

type,  B-B  and  C-C  (V2  a/2),  is  greater  than  that  of  the  B-C  pair  (a/2).  Then,  when  the 
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system  relaxes,  two  solutes  of  different  type  will  have  more  chance  to  enter  into  contact 
that  two  solutes  of  the  same  type  as  the  composition  increases.  This  is  reflected  in  the 
behaviour  of  Z(  li,  which  starts  to  increase  from  zero  at  X^=X(  ~0. 1 0,  while  Zbb  and  Zee 
start  at  Xb=Xq~0.  1 5 . 

If  we  require  that  the  optimal  composition  is  again  that  which  retains  chemical 
ordering,  then  the  maximum  density  (note  that  both  solutes  are  smaller  than  the  solvent 
and  so  the  density  increases  with  solute  composition)  for  which  all  types  of  solute-solute 
contact  are  avoided  is  roughly  0.08  for  each  solute,  or  the  total  solute  composition  to  be 
0.16,  close  to  the  value  we  found  for  the  binary  system  with  a  small  solute. 


(a) 
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(b) 


Figure  6:  Coordination  numbers  Zy  for  the  ternary  system  studied,  RpJR^O.HO  and  Rc/Ra=0.62:  a)  Zaa,  ZBA 
and  ZCA  and  b)  ZBB,  ZCB  and  Zcc,  the  inset  shows  in  more  detail  the  behaviour  in  the  low  composition 
region.  The  full  line  and  symbols  (dashed  line)  correspond  to  the  constrained  (unconstrained)  system. 


2)  Comparison  of  the  MPA  with  a  chemically  ordered  amorphous  packing 

The  most  direct  test  of  the  MPA  is  to  compare  it  with  the  analogous  amorphous  structure 
with  chemical  order.  We  have  generated  such  structures  by  starting  the  solutes  in  a 
random  close  packed  structure  instead  of  on  fee  lattice  sites.  As  described  in  the 
following  section,  this  algorithm  produces  chemically  ordered  structures  for  XB  <0.18. 
While  this  comparison  has  been  complicated  by  the  appearance  of  the  fully  crystalline 
states  at  some  size  ratios,  we  have  obtained  some  useful  results.  In  Table  2  we  compare 
the  value  of  ZBa  for  XB  <  0.2  for  the  MPA  and  amorphous  systems 


Pr/Pa 

Zba(MPA) 

ZBA( amorphous) 

0.62 

8 

7.9 

0.80 

10 

9.8 

Table  2.  A  comparison  of  the  values  of  ZBA  found  in  the  MPA  and  a  chemically  ordered  amorphous 
packing  of  a  binary  mixture  of  sphere  with  the  values  of  RB/RA  indicated. 

These  preliminary  results  indicate  that  the  MPA  is  successfully  describing  important 
features  of  the  local  structure  of  the  amorphous  state. 


3)  Establishing  the  compositional  limit  on  chemical  order  in  amorphous  packings  of 
mixtures  of  spheres. 
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Starting  from  either  the  MPA  or  an  amorphous  configuration,  we  have  identified  a 
solute  composition  above  which  chemical  order  breaks  down  as  signaled  by  the  increase 
in  ZBb.  The  existence  of  this  limiting  composition  is  of  importance  in  predicting  the 
optimal  composition  for  glass  formation  in  an  alloy.  As  part  of  this  project  we  have 
formulated  the  following  proposition: 

The  optimal  composition  for  glass  stability  is  that  which  produces  the  maximum 
density  while  a)  not  allowing  for  solute-solute  contact,  and  b)  not  resulting  in  isolated 
solvents,  i.e.  solvents  without  a  solute  neighbour. 

As  shown  in  Figure  6,  criteria  a)  and  b)  from  the  stability  proposition  coincide  at  a 
fairly  well  defined  composition,  a  composition  that  depends  on  the  choice  of  radius  ratio. 


(a)  (b) 

Figure  6:  fraction  of  ‘isolated  A’  spheres  not  having  contact  with  any  B  and  the  fraction  of  ‘B-B  contacts’ 
as  a  function  of  XB  for  both  system  size  ratios:  a)  RB/RA=0.80  and  b)  -/?b/^a=1.12 


We  have  carried  out  a  systematic  study  of  the  composition  dependence  of  Zbb  in  the 
amorphous  packings,  the  results  of  which  are  shown  in  Figure  7a.  In  Figure  7b  we  plot 
the  observed  dependence  of  this  upper  compositional  limit  c*  on  the  size  ratio.  We  find  a 
relatively  weak  dependence  of  c  on  the  size  ratio  with  0.14  <  c  <  0.19  over 
0.6  <  Rb  /  R4  <  1 .25  .  c  decreases  steadily  with  increasing  size  of  the  solvent  since  it  takes 
an  increasing  number  of  solvent  particles  per  solute  to  keep  the  solutes  apart. 
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Figure  7:  a)  Solute-solute  coordination  number  ZBB  as  a  function  of  XB  for  each  size  ratio  Rb/Ra,  b) 
estimated  upper  limit  for  XB  as  a  function  of  Rb/Ra  (chemically  ordered  systems  with  the  B  spheres  in 
random  close  packed  positions)  corresponding  to  a).  All  the  lines  are  guides  to  the  eye. 


4)  Comparison  of  chemically  ordered  amorphous  sphere  packing  with  a  model 
metal-metalloid  mixture 

It  is  of  interest  to  compare  the  coordination  and  structure  of  the  constrained  hard  sphere 
mixture  and  a  more  realistic  model  with  analogous  size  ratio.  Here  we  consider  a  binary 
mixture  of  particles  interacting  by  a  Lennard-Jones  (LJ)  potential  of  the  fonn 
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where  lengths  cAa  =1.0,  oAB  =  0.8  and  gBb=0.88.  The  interaction  strengths  are  sAa  =1.0, 
sAb=  1.5  and  Sbb  =  0.5  so  we  expect  strong  chemical  ordering.  This  model  has  been 
studied  extensively  in  the  context  of  the  glass  transition  and,  in  that  field,  is  referred  to  as 
the  KA  mixture.  Considering  just  the  AA  and  AB  lengths,  this  corresponds  to  a  hard 
sphere  size  ratio  of  0.6.  A  better  approximation  would  involve  the  positions  of  the  first 
peaks  in  gAA  and  gAB:  for  the  quenched  KA  mixture  it  is  obtained  an  effective  size  ratio 
of  RB/RA~0.65  (rAAmax~1.07oAA  and  rABmax~0.88  oAA).  Then,  we  can  compare  the 
quenched  KA  and  the  hard-sphere  system  with  RB/RA=0.62  (corresponding  to  a 
theoretical  ZBA=8).  We  use  the  relationship  RA+RB~r ABmax  to  scale  distances  from  Ra 
units  into  oAA.  Again,  this  scaling  is  possible  because  there  is  almost  no  B-B  contacts  in 
the  KA  system. 
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Figure  8:  Average  coordination  numbers  ZHA  and  ZBB  for  quenched  states  of  the  KA  system  (T=0.4)  and 
two  hard  sphere  mixtures  with  short  range  order  (RB/R a=0.62  and  0.80)  as  a  function  of  the  composition 
XB. 


In  Figure  8  we  plot  the  coordination  numbers  ZAb  and  ZBb-  The  chemically  ordered  hard 
sphere  mixture  does  a  quite  reasonable  job  of  predicting  the  coordination  numbers  and 
their  concentration  dependence  in  the  LJ  model.  We  find  the  LJ  results  lie  a  little  above 
the  hard  sphere  Zab  and  Zbb-  BB  contacts,  explicitly  excluded  at  low  concentration  in  the 
hard  sphere  mixture,  are  present  in  small  numbers  in  the  unconstrained  LJ  system. 

A  sterner  comparison  of  the  two  structures  is  provided  by  the  radial  distribution 
functions.  In  Figure  9  we  compare  the  solute-solute  distribution  gBB(r)  for  the  hard  sphere 
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and  LJ  mixture.  The  LJ  mixture  shows  considerable  more  well  defined  structure  in  the 
form  of  peaks  at  a  number  of  distances.  The  HS  mixture,  in  contrast,  has  a  single  peak  (in 
a  position  similar  to  a  peak  in  the  LJ  mixture)  with  broad  shoulders.  As  demonstrated 
previously,  the  individual  peaks  in  the  LJ  mixture  can  be  assigned,  unequivocally,  to 
contributions  from  different  types  of  BAm  bonds.  We  find  that  the  major  peak  in  the  HS 
mixture  and  the  peak  at  the  similarly  long  distance  in  the  LJ  mixture  arises  from  the  BA1 
bonds  (a  shared  vertex  of  the  solute  coordination  shells).  We  see  that  missing  structure  in 
the  HS  mixture  is  due  to  the  almost  complete  absence  of  well  defined  triangular  face 
sharing  (BA3)  and  a  very  broad  peak  associated  with  the  BA2  (edge  sharing)  bond.  The 
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Figure  9:  Partial  (full  line)  and  total  (dotted  line)  solute-solute  correlation  function  for  hard 
spheres  (HS)  with  RB/RA= 0.62  and  quenched  KA  system  from  T=0.4,  both  at  Ab=0.20:  a) 
BA1,  b)  BA2,  c)  BA3  and  d)  BA4.  The  distance  axis  for  the  HS  system  has  been  scaled 
down  by  a  factor  0.54  in  order  to  change  units  from  RA  to  Oaa- 
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impression  is  that  the  coordination  shells  of  the  solute  in  the  HS  mixture  are  considerably 
less  organised  than  the  LJ  mixture.  While  our  constraint  has  provided  chemical  order,  we 
cannot  duplicate  the  strong  local  correlations  arising  from  the  particle  interactions. 
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